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Abstract 

The quark mass dependence of nucleon magnetic moments is explored over 
a wide range. Quark masses currently accessible to lattice QCD, which lie 
beyond the regime of chiral perturbation theory (%PT), are accessed via the 
cloudy bag model (CBM). The latter reproduces the leading nonanalytic be- 
havior of xPT, while modeling the internal structure of the hadron under in- 
vestigation. We find that the predictions of the CBM are succinctly described 
by the simple formula, //2v(m,r) = £% /(I + ctm^ + (3m%), which reproduces 
the lattice data, as well as the leading nonanalytic behavior of xPT. As this 
form also incorporates the anticipated Dirac moment behavior in the limit 
m n — > oo, it constitutes a powerful method for extrapolating lattice results 
to the physical mass regime. 
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I. INTRODUCTION 



The fundamental theory of the strong interactions, Quantum Chromodynamics (QCD), 
is a nonperturbative field-theory in which the self-coupling of the gauge field leads to a non- 
trivial vacuum structure. The complexity of the QCD vacuum is manifest in nonvanishing 
vacuum expectation values of quark and gluon operator products having vacuum quantum 
numbers. The only known way to directly calculate the properties of QCD is through the 
formulation of a lattice gauge theory where the fields are described on a discrete space-time 
lattice. 

The lattice formulation of QCD is well established Recent advances in lattice action 
improvement on anisotropic lattices are greatly facilitating the reduction of systematic un- 
certainties associated with the finite lattice volume and the finite lattice spacing. However, 
direct simulation of QCD for light current quark masses, near the chiral limit, remains com- 
putationally intensive. As such, the present approach of calculating the properties of QCD 
using quark masses away from the chiral regime and extrapolating to the physical world is 
likely to persist for the foreseeable future. 

The difficulty associated with this approach is illustrated by the rapid rise of the pseu- 
doscalar mass for small increases in the quark mass away from the chiral limit as governed 
by 

m -- 2ma «, (!) 

J n 

where q denotes a light quark. Typical quark masses considered in lattice simulations place 
~ 500 MeV, a value significantly larger than the physical pion mass, fi = 139.6 MeV. A 
mass the order of 500 MeV is sufficient to largely suppress the pion cloud contribution to 
hadronic observables, whereas near the chiral limit, the pion cloud can make a significant 
contribution to them [0] . Hence it is imperative to extrapolate any observable simulated on 
the lattice to the physical world using a function motivated by the physics of the pion cloud. 
Since the pion cloud contributions are small in lattice simulations, one has to look beyond 
the lattice simulation results at present. 

Historically, lattice results were often linearly extrapolated with respect to m\ to the 
chiral limit, particularly in exploratory calculations. More recently the focus has turned to 
chiral perturbation theory (xPT), which provides predictions for the leading nonanalytic 
quark- mass dependence of observables in terms of phenomenological parameters pHTTf. 
Application of to the extrapolation of lattice simulation data is now standard for 



hadron masses and decay constants |I] . However, earlier attempts (12| to apply %PT predic- 
tions for the quark-mass dependence of baryon magnetic moments failed, as the higher order 
terms of the chiral expansion quickly dominate the truncated expansion as one moves away 
from the chiral limit. To one meson loop, expresses the nucleon magnetic moments as 

EMI 

/J-n = Ak) + Cim T + c 2 ml log ml + c 3 m 2 n H , (2) 

where /io and C3 are fitted phenomenologically and c\ and C2 are predicted by xPT. The 
m\ log m\ term quickly dominates as m n moves away from the chiral limit making contact 
with the lattice results untenable. 
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Lattice QCD results for baryon magnetic moments fT5}-|T7|| remain predominantly based 
on linear quark mass (or m 2 ) extrapolations of the moments expressed in natural magnetons. 
This approach systematically underestimates the measured moments by 10 to 20%. Finite 
lattice volume and spacing errors are expected to be some source of systematic error. How- 
ever, xPT clearly indicates the linear extrapolation of the simulation results is also suspect. 
As such, it is imperative to find a method which can bridge the void between the realm of 
%PT and lattice simulations. 

We report such a method, which provides predictions for the quark mass (or m 2 ) depen- 
dence of nucleon magnetic moments well beyond the chiral limit. In particular, we use the 
cloudy bag model (CBM), which involves a relativistic quark model (the MIT bag) coupled 
to the pion field in such a way as to restore chiral symmetry |18|,|19| . The corresponding pion 
loop corrections to physical observables reproduce the leading non-analytic behavior of %PT. 
However, the loop which gives rise to the leading non-analytic behavior involves two pion 
propagators. Because the loop is regulated by a form factor, related to the finite size of the 
hadron, its contribution is suppressed like l/m^ as becomes large. This feature makes it 
possible to address the larger quark masses simulated in lattice QCD in a convergent way. 

With some tuning of the bag radius, the form factor at the nNN vertex and the current 
quark mass, within the framework of the CBM, we find that it is possible to reproduce 
the lattice QCD simulations of nucleon magnetic moments and smoothly extrapolate to 
the experimentally measured values. We then propose a simple phenomenological relation 
designed to reproduce the leading nonanalytic structure of ;^PT and provide the Dirac- 
moment mass dependence in the heavy quark-mass regime. Finally, we illustrate how such 
a relation can be used in future lattice QCD calculations. These results may also be useful 
in clarifying issues surrounding the higher order terms of the chiral expansion of x?T. 



II. LATTICE QCD SIMULATION DATA 

We consider two independent lattice simulations of the nucleon electromagnetic form 
factors. Both calculations employ three-point function based techniques [^0| utilizing the 



conserved vector current such that no renormalization is required in relating the lattice 
results to the continuum. Ref. utilizes twenty-eight quenched gauge configurations on a 
24 x 12 x 12 x 24 periodic lattice at (3 = 5.9, corresponding to a lattice spacing of 0.128(10) 
fm. Moments are obtained from the form factors at 0.16 GeV 2 by assuming equivalent q 2 
dependencies for the electric and magnetic form factors. Ref. |T(| utilizes twelve quenched 
gauge configurations on a 16 3 x 24 periodic lattice at /3 — 6.0. The nucleon mass |21| 



corresponds to a lattice spacing of 0.091(3) fm. Moments are obtained from dipole fits 
to the form factors. Uncertainties are statistical in origin and are estimated by a single 



elimination jackknife fl22| . Despite having different lattice volumes and lattice spacings, the 
results from the two calculations agree well within errors. However, experience suggests 
that the lattice spacings and volumes used in these investigations may give rise to scaling 
violations and finite size effects the order of 15 to 20% from the infinite volume continuum 
limit. 

Ideally, one would like to perform the analysis of chiral corrections, using results from 
full QCD, as opposed to quenched QCD. Unfortunately, such results are not yet available. 
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Instead we will utilize these results under the standard approach of correcting the lattice 
scale by fixing the lattice spacing using the nucleon mass. In the absence of any known way 
to correct for the nonperturbative effects of quenching, we will assume that these lattice 
results are a reasonable representation of full QCD. Since the quark masses simulated on 
the lattice are somewhat heavy, the dominant contributions from quark loops neglected in the 
quenched approximation are largely perturbative and accounted for in the renormalization 
of the lattice spacing. As we shall see, the errors due to quenching are likely the same order 
of magnitude as the statistical uncertainties. 



III. THE CLOUDY BAG MODEL 



The linearized CBM Lagrangian with the pseudoscalar pion-quark coupling (to order 
1/U ) is given by ppl 



£ = [qiirfdp - m q )q - B}9 V - -yyos 

+ \( d ^) 2 ~ \ m W ~ T^qisT ■ izq5 s , (3) 

where B is the bag constant, is the 7r decay constant, 9y is a step function (unity inside 
the bag volume and vanishing outside) and 5s is a surface delta function. In a lowest order 
perturbative treatment of the pion field, the quark wave function is not affected by the pion 
field and is simply given by the MIT bag solution P5|| . Our calculation is carried out in the 
Breit frame with the center-of-mass correction for the bag performed via Peierls-Thouless 
projection. The detailed formulas for calculating nucleon electromagnetic form factors in 
the CBM are given in Ref. |4|. 



In the CBM, a baryon is viewed as a superposition of a bare quark core and bag plus 
meson states. Both the quark core and the meson cloud contribute to the baryon magnetic 



moments. These two sources are balanced around a bag radius, R = 0.7 — 1.1 fm [25]. A 
large bag radius suppresses the contributions from the pion cloud, and enhances the contri- 
bution from the quark core. The minimal coupling principle is used for the electromagnetic 
interaction. The nucleon magnetic moments can be written as //jy = Gm(0) + G^(0), where 
Gm is due to ^yqq coupling and from 'jtttt coupling. To one loop, the CBM reproduces 
the leading non-analytic behavior of x?T. The processes included in this calculation are 
illustrated in Fig. [[]. 

For the ttNN vertex, we take the conventional ttNN coupling constant, f% NN = 0.0771. 
Instead of the generic form, the nNN form factor is replaced by the phenomenological, 
monopole form, u(k) = (A 2 — yU 2 )/(A 2 + k 2 ), where k is the loop momentum and A is a 
cut-off parameter. 

In the standard CBM treatment, where the pion is treated as an elementary field, the 
current quark mass, m q , is not directly linked to m n . Most observables are not sensitive to 
this parameter, as long as it is in the range of typical current quark masses. For our present 
purpose it is vital to relate the m q inside the bag with m n . Current lattice simulations 



indicate that m 2 is approximately proportional to m q over a wide range of quark masses [p26 



Hence, in order to model the lattice results, we scale the mass of the quark confined in the 
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FIG. 1. Schematic illustration of the processes included in the CBM calculation. B and C 
denote intermediate state baryons and include N and A. 

bag as m q = (m 7r ///) 2 TOg , with mf being the current quark mass corresponding to the 
physical pion mass, mf is treated as an input parameter to be tuned and lies in the range 
6 to 7 MeV. 

The parameters of the CBM are obtained as follows. The bag constant B and the 
phenomenological parameter zq are fixed by the physical nucleon mass and the stability 
condition, dM^/dR = 0, for a given Rq and mf\ For each subsequent value of the pion 
mass or the quark mass considered, u and R are determined simultaneously from the 
linear boundary solution of the bag EBJ and the stability condition, which provides R A = 



(3uo — Zo) I (AttB). Using the lattice data and the experimental measurement, the parameters 
Ro, A and mf are tuned to reproduce the experimental moment while accommodating the 
lattice data. 

Fig. shows the mass dependence of the bag model properties and the pion-cloud con- 
tribution G\j{Q). The bare bag probability, Z 2) the quark ground state frequency, u , and 
the bag radius, R, are plotted as a function of m^. As m n increases, the pion-cloud contri- 
bution decreases very quickly and becomes quite small for large quark masses - especially 
in the range corresponding to the current lattice calculations. On the other hand, the bag 
properties evolve relatively slowly. As a result, the dominant influences of Rq, A governing 
the irNN coupling, and rriq are located in separate regions of m n . The magnitude of the 
magnetic moments in the small region is controlled by A and Rq, while a variation of 
rriq is more effective in the large m n region, where the pion cloud nearly vanishes. 

IV. DISCUSSION OF RESULTS 

The nucleon magnetic moments calculated in the CBM are shown in Fig. for the proton, 
and Fig. f|for the neutron by the solid line. The lattice results (• |[L5|| , ■ [H)]) are also plotted 



in these figures. It is possible to simultaneously reproduce the existing lattice simulation 
results for the nucleon magnetic moments at large m n as well as their physical values using 
CBM parameters within previously established ranges. These parameters are summarized 
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FIG. 2. The pion contribution to the proton magnetic moment G^(0), and the properties of 
the bag including the bare bag probability Z2, the ground state frequency ujq, and the bag radius 
R in variation with the pion mass. 

in Table |. The dashed lines in Figs. |3| and |] indicate the corresponding results for the MIT 
bag model. There, without the pion cloud, the m 2 dependence of the magnetic moments 
becomes nearly linear. This clearly shows the significance of the meson cloud, especially in 
the small m n regime. 

Since the pion cloud in the quenched approximation is quite different from that of full 
QCD PH5|.f27|| , one might regard the difference between the CBM and MIT model results 
(i.e. the pion contribution) as indicative of the absolute systematic error associated with the 
quenched approximation. For the quark masses actually simulated on the lattice, this error 
is the order of 15%. However, some of this difference is already accounted for in the lattice 
results through a renormalization of the lattice spacing used in expressing the moments in 



TABLE I. CBM parameters and optimal fit parameters for the fit function of Eq. (|j). Other 
CBM parameters are Zq = 2.59 and B 1 ^ = 144 MeV. At physical pion mass, the experimental 
magnetic moments (2.79 and -1.91 for the proton and neutron) are reproduced in the CBM. 



N 


R (fm) 


A (GeV) 


m { g ] (MeV) 


(0) 

/% 


a (GeV" 1 ) 


(GeV" 2 ) 


proton 


1.0 


0.68 


4.8 


3.31 


1.37 


0.452 


neutron 


1.0 


0.59 


4.8 


-2.39 


1.85 


0.271 
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FIG. 3. The proton magnetic moment as calculated in lattice QCD (• LDW Ref. [15], ■ WDL 
Ref. fl6|]), the cloudy bag model (CBM) which includes the pion cloud contribution and the MIT 
bag model (MIT) where the pion cloud contribution of the CBM is omitted. Also illustrated is a 
fit of the simple analytic form of Eq. (0) to the CBM results. The point at the physical value of 

is the experimental measurement and is used to constrain the parameters of the CBM. 



nuclear magnetons. Hence, for the purposes of this investigation, it is reasonable to accept 
the quenched results as an approximate representation of the full QCD result. 

V. ENCAPSULATING FORM 



Having established the quark mass dependence of the nucleon moments over a very wide 
range, we now turn to encapsulating these results in a simple analytic form that might be 
used in future lattice QCD extrapolations of simulation results. The following function 



(o) 



1 + am, + (3m\ 



(4) 



is matched to the CBM model results by tuning the three parameters fiff , a and (3. These 
fits are also illustrated in figures ^ and |4j. This functional form provides the correct limiting 
behavior as a function of m n . As m v — > 0, Eq. (|j) may be expanded as 



[i - 



am. 



[a 



0)< + ■■■} 



(5) 



such that the leading nonanalytic behavior is proportional to as required by xPT. For 
large m n , Eq. (B) leads to 



7 



0.0 




-2.0 1 1 1 1 1 1 1 1 ' 

0.0 0.5 1.0 1.5 2.0 

m, 2 (GeV 2 ) 



FIG. 4. The neutron magnetic moment. Symbols and lines are as described for Fig. ||. 



Mm,) = -^(l-^ + --) , (6) 

such that the magnetic moments decrease as l/m q for increasing quark mass, precisely as the 
Dirac moment requires. In short, this function reproduces the leading nonanalytic behavior 
of xPT, provides the desired mass dependence of the Dirac moment in the heavy quark mass 
limit and also has the required shape in our region of interest. The resulting parameters for 
this fit are also listed in Table |. 

Because the leading nonanalytic terms predicted in xPT for baryon magnetic moments 
are proportional to the pseudoscalar meson mass, the kaon cloud is often regarded as the most 
important and dominant contribution. Instead, in the CBM with its natural high momentum 
cut-off, the overall contribution of the kaon loop is strongly suppressed. Moreover, the kaon 
cloud contribution to the nucleon magnetic moments will not display the same significant 
curvature associated with the pion cloud. In consequence, the kaon cloud contributions may 
simply be absorbed into the fit parameters fi^ and (3. 

To evaluate the ability of Eq. (^) to encompass the CBM predictions, the lattice QCD 
data and maintain the correct leading nonanalytic behavior of xPT, we compare the coef- 
ficient of m n in Eq. (|5J) with that predicted by xPT. Table Q summarizes numerical values 
for this coefficient. The one-loop corrected estimates for the D and F coefficients provide 
better agreement between %PT and experiment for many observables [I3|. The similarity 



between this xPT estimate and the coefficient from the encapsulating form is encouraging. 
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TABLE II. Values for the coefficient of the leading nonanalytic term obtained from the chiral 
expansion of Eq. (Q) adjusted to fit the CBM results for the proton and neutron and from xPT. 
H$ and a are from Table |. Dq and Fq are tree-level coefficients (Dq + Fq = gA = 1.27) while D\ 
and F\ are the one- loop corrected estimates of Ref. [13|. Units are GeV -1 . In the column headings, 
the upper and lower signs correspond to the proton and neutron respectively. 



Nucleon 


Encapsulating Form 


Chiral Perturbation 


Theory 




(0) 


m N {D + F ) 2 


m N {D 1 + F l ) 2 


proton 


-4.54 


-6.97 


-4.41 


neutron 


4.42 


6.97 


4.41 



VI. EXTRAPOLATION OF LATTICE QCD DATA 



Future lattice QCD studies of octet baryon magnetic moments will make better con- 
tact with experiment by adopting the one-loop corrected coefficient of from %PT and 
performing a two parameter fit of Eq. (H) to the simulation data [p8[| . The utility of this 



approach is illustrated in Fig. g. The nucleon magnetic moments at the physical pion mass 
obtained from this extrapolation are 

Hp = 2.85(22) fi N and pi n = -1.90(15) fi N (7) 

and agree surprisingly well with the experimental measurements, 2.793 and —1.913 
respectively. The fit parameters (fx^°\ P) in units of (hn, GeV -2 ) are (3.39(23), 0.58(16)) 
and (—2.40(16), 0.41(16)) for the proton and neutron respectively. We note that the data 
required to do a fit of the lattice results in which covariances are taken into account is no 
longer available. As such, the uncertainties quoted here should be regarded as indicative 
only. 



VII. CONCLUSIONS 



In summary, we have explored the quark mass dependence of nucleon magnetic moments. 
Quark masses beyond the regime of chiral perturbation theory have been accessed via the 
cloudy bag model which reproduces the leading nonanalytic behavior of x?T and provides 
internal structure for the hadron under investigation. We find that the predictions of the 
CBM are succinctly described by a simple formula which reproduces the leading nonanalytic 
behavior of xPT in the limit m % —>■ and provides the anticipated Dirac moment behavior in 
the limit m n — ► oo. The significance of nonlinear behavior in extrapolating nucleon magnetic 
moments as a function of m q to the chiral regime has been evaluated. We find that the leading 
nonanalytic term of the chiral expansion dominates from the chiral limit up to the physical 
pion mass. Beyond the physical mass, higher order terms become important and dominate. 
This curvature, neglected in previous linear extrapolations of the lattice data, can easily 
account for the departures of earlier lattice estimates from experimental measurements. As 
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FIG. 5. Extrapolation of lattice QCD magnetic moments (• LDW Ref. ||, ■ WDL Ref. @) 
for the proton (upper) and neutron (lower) to the chiral limit. The curves illustrate a two parameter 
fit of Eq. (|j) to the simulation data in which the one-loop corrected chiral coefficient of is taken 
from xPT. The experimentally measured moments are indicated by asterisks. 

finite volume and lattice spacing artifacts are eliminated in future lattice QCD simulations, 
it will be interesting to see if the fit parameters adjust accordingly to maintain and perhaps 
improve the level of agreement seen in this investigation. We advocate the use of the function 
in Eq. (Q) in future lattice QCD investigations of octet baryon magnetic moments. 
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